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A particular class of space-time, with a tachyon field, (f}, and a barotropic fluid constituting the 
matter content, is considered herein as a model for gravitational coUapse. For simplicity, the tachyon 
potential is assumed to be of inverse square form i.e., V{4>) ~ 4'~'^- Our purpose, by making use of the 
specific kinematical features of the tachyon, which are rather different from a standard scalar field, is 
to establish the several types of asymptotic behavior that our matter content induces. Employing a 
dynamical system analysis, complemented by a thorough numerical study, we find classical solutions 
corresponding to a naked singularity or a black hole formation. In particular, there is a subset where 
the fluid and tachyon participate in an interesting tracking behaviour, depending sensitively on the 
initial conditions for the energy densities of the tachyon fleld and barotropic fluid. Two other classes 
of solutions are present, corresponding respectively, to either a tachyon or a barotropic fiuid regime. 
Which of these emerges as dominant, will depend on the choice of the barotropic parameter, 7. 
Furthermore, these collapsing scenarios both have as final state the formation of a black hole. 

PACS numbers: 04.20.Dw, 04.70.Bw, 04.20.Cv, 97.60.Lf 



I. INTRODUCTION 



The study of the final state of the gravitational collapse of initially regular distributions of matter is one of the open 
problems in classical general relativity, having attracted remarkable attention in past decades. When a sufficiently 
massive star exhausts all the thermonuclear sources of its energy, it would undergo a collapsing scenario due to its own 
gravity, without reaching a final state in terms of a neutron star or white dwarf. Under a variety of circumstances, 
singularities will inevitably emerge (geodesic incompleteness in space-time), matter densities and space-time curvatures 
diverging. Albeit the singularity theorems [l][2] state that there exist space-time singularities in a generic gravitational 
collapse, they provide no information on the nature of singularities: the problem of whether these regions are hidden 
by a space-time event horizon or can actually be observed, remains unsolved. The cosmic censorship conjecture 
(CCC), as hypothesized by Penrose [3], conveys that the singularities appearing at the collapse final outcome must 
be hidden within an event horizon and thus no distant observer could detect them. A black hole forms. Although 
the CCC plays a crucial role in the physics of black holes, there is yet no proof of it, due to the lack of adequate 
tools to treat the global characteristics of the field equations. Nevertheless, in the past 30 years many solutions to the 
field equations have been discovered, which exhibit the occurrence of naked singularities, where the matter content 
has included perfect and imperfect fiuids [4j [5] , scalar fields [6] , self-similar models [7] and null strange quarks [8] . 
Basically, it is the geometry of trapped surfaces that decides the visibility or otherwise of the space-time singularity. 
In case the collapse terminates into a naked singularity, the trapped surfaces do not emerge early enough, allowing 
(otherwise hidden) regions to be visible to the distant observers. 

The gravitational collapse of scalar fields is of relevance [9, , owing to the fact that they are able to mimic other 
types of behaviours, depending on the choice of the potentials. Scalar field models have been extensively examined 
for studying CCC in spherically symmetric models [lOj . non-spherically symmetric models |llj and also for static 
cases [H] . Their role in understanding the machinery governing the causal structure of space-time was available since 
the 90's, when the numerical solutions exhibiting naked singularities were found numerically by Choptuik [13^ and 
analytically by Christodoulou il4j . There are in the literature a few papers discussing gravitational collapse in the 
presence of a scalar field joined by a fluid for the matter content [T5l [TB]: in summary, a black hole forms in these 
collapsing situations. However, to our knowledge, a tachyon scalar field has not yet been considered regarding whether 
a black hole or naked singularity forms, that is to say, in the CCC context, together with a fiuid. Tachyon fields arise in 
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the framework of string theory [T7] and have been of recent use in cosmology [IB] ■ The action for the tachyon field has 
a non-standard kinetic term |19j , enabling for several effects whose dynamical consequences are different from those of 
a standard scalar field [SO]- Namely, other (anti-)friction features that can alter the outcome of a collapsing scenario. 
This constitutes a worthy motivation to investigate the scenario where a tachyon field is added to a barotropic fluid, 
both constituting the matter content present in the collapse process: on the one hand, the fluid will play the role of 
conventional matter from which a collapse can proceed into, whereas, on the other hand, the tachyon would convey, 
albeit by means of a simple framework, some intrinsic features from a string theory setting. Restricting ourselves 
herein to the tachyon as the intrinsic string ingredient influencing the collapse, let us nevertheless point that many 
other string features could be incorporated in subsequent similar studies PjUZlIIl]. Our purpose, in this paper, by 
investigating the gravitational collapse of a barotropic fluid together with a tachyon field, is therefore to establish 
the types of final state that can occur (i.e., whether a black hole or a naked singularity emerges, in the context of 
the CCC), which matter component will determine the outcome. In particular, if the late time tachyon behaviour, 
possibly competing with the fluid and eventually becoming dominant, could allow interesting features to appear. 

We then organize this paper as follows. In section|ll]we give a brief revie w on the gravitational collapse of a specific 
space-time, namely the marginally bounded case (cf. \'21\). In section III we study, by means of a dynamical system 
analysis, the gravitational collapse employing a tachyon and a barotropic fluid as the matter content. The analytical 
study is complemented by a careful numerical investigation. In section IV we present our conclusions and a discussion 
of our results. 



II. A BRIEFING ON MARGINALLY BOUND GRAVITATIONAL COLLAPSE 

In this section, we will discuss the space-time region inside the collapsing sphere which will contain the chosen matter 
content. An isotropic Friedmann-Robertson- Walker (FRW) metric, in comoving coordinates, will be considered as the 
interior space-time for the gravitational collapse. However, in order to study the whole space-time, we must match 
this interior region to a suitable exterior. In the model herein, it is convenient to consider a spherically symmetric and 
inhomogeneous space-time such as the Schwarzschild or the generalized Vaidya geometries to model the space-time 
outside the collapsing sphere. 

The interior space-time for the marginally bounded case (cf. ^21j) can be parametrized as 

ds^ = -dt^ + a(tfdr^ + R{t, rfdn'^, (2.1) 

where t is the proper time for a falling observer whose geodesic trajectories are labeled by the comoving radial 
coordinate r, R{t, r) = ra[t) is the area radius of the collapsing volume and dVl^ being the standard line element on 
the unit two-sphere. The Einstein's field equations for this model can be presented [I] as (we have set SttC? = c = 1 
throughout this paper) 

F —F 

(2.2) 



^ Rm^r' " R^R 
where the mass function 

F{t,r) = RR^, (2.3) 

is the total gravitational mass within the shell, labeled by the comoving radial coordinate r with ", r = dr" and 
"• = 9(" . Since we are interested in a continuous collapsing scenario, we take a < 0, implying that the area radius of 
the collapsing shell, for a constant value of r, decreases monotonically. Splitting the above metric for a two-sphere 
and a two dimensional hypersurface, normal to the two-sphere, we have 

ds^ haidx^-dx" + R{t, rfdfl^, hat = diag [-1, a{t)^] , (2.4) 

whereby the Misner-Sharp energy [32] is defined to be 
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Eit, r) = [1 _ h'^^d^Rd.R] ^ ^ . (2.5) 



To discuss the final state of gravitational collapse it is important to study the conditions under which the trapped 
surfaces form. From the above definition, it is the ratio 2E{t, r)/R{t, r) that controls the formation of trapped surfaces 
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[25] (note that the mass function defined above is nothing but twice the Misner-Sharp mass). Introducing the null 
coordinates 



1 



1 



= -= [dt - a{t)dr] , = - ^ [dt + a{t)dr] , 



V2 



V2 



the metric (2.4) can be recast into the double-null form 

ds^ = -2d^+dC +R{t,rfdn'^. 



(2.6) 



(2.7) 



The radial null geodesies are given by ds^ = 0. We then may deduce that there exists two kinds of null geodesies 
which correspond to = constant and = constant, the expansions of which are given by 



0± = ^d±R, 



where 



d 



-V2 



dt 



dr 



d 

of 



-V2 



dt 



a{t) 



(2., 



(2.9) 



The expansion parameter measures whether the bundle of null rays normal to the sphere is diverging {Q± > 0) or 
converging (0± < 0). The space-time is referred to as trapped, untrapped and marginally trapped if 



e+e_>o, e+e_<o, e+e_ = o, 



(2.10) 



respectively, where we have noted that h''''daRdbR = %+^' . The third case characterizes the outermost boundary 
of the trapped region, the apparent horizon equation being R^ — 1. From the point of view of equation (2.5), we see 
that at the boundary of the trapped region, 2E{t, r) = R{t, r). Thus, in the regions where the mass function satisfies 
the relation F{t,r) > R(t,r), the trapping of light does occur, whereas for the regions in which the mass function is 
less than the area radius i.e., F(t,r) < R{t,r), no trapping happens [22j . In other words, in the former there exists 
a congruence of future directed null geodesies emerging from a past singularity and reaching to distant observers, 
causing the singularity to be exposed (a naked singularity forms). In the latter, no such families exist and trapped 
surfaces will form early enough before the singularity formation: regions will be hidden behind the event horizon, 
which lead to the formation of a black hole [T] (a similar interpretation has been given in PSI 124) '). 

In order to further illustrate this specific gravitational collapse process, let us employ a very particular class of 
solutions. We assume that the behavior of the matter density, for a <C 1, is given by the following ansatz^ (cf. ref. 



Pit) 



(2.11) 



where c and n are positive constants^. Integration of the first equation in (2.2) gives the following relation for the 
mass function as 



(2.13) 



Hence we have 



2 2~n 



R = r'' 



(2.14) 



It is therefore possible to identify the following outcomes: 



^ Notice that, since the interior space-time |2.4| is spatially homogeneous, it is obvious that, the energy-momentum tensor must also be 

spatially homogeneous. 

^ Using the energy density ( |2.1l[ | and pressure l |2.15| l, the weak energy condition can be written as follows 

P + P=^>0- (2-12) 

The weak energy condition is satisfied throughout the collapse process (see also next section for more details regarding the issue of the 
energy conditions). 
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Figure 1: The potential V = Vo(f> ^, denoting tlie </> < and </> > branclies as well as the asymptotic stages. 



• For < ri < 2, the ratio F/R < 1 is less than unity throughout the collapse and if no trapped surfaces exist 
initially, due to the regularity condition {F{tQ, r)/R{tQ, r) < 1 for a suitable value of c and r), then no trapped 
surfaces would form until the epoch a{ts) — 0. More precisely, there exists a family of radial null trajectories 
emerging from a naked singularity. 

• For n > 2, the ratio F/R goes to infinity and trapped surfaces will form as the collapse evolves, which means 
that singularity will be covered and no radial geodesies can emerge from it. Thus, a black hole forms. 



From equation (2.11 ), we may easily find a relation between the energy density and pressure as 

ri — 3 



(2.15) 



which shows that for n < 2, where a naked singularity occurs, the pressure gets negative values f27'. From equations 
(2.3), (2.11 1 and bearing in mind the continuous collapse process (a < 0) we can solve for the scale factor as 



a{t) 



(to - t) 



(2.16) 



where to is the time at which the energy density begins increasing as a " and Oq corresponds to tQ. Assuming that 
the collapse process initiates at to — 0, the time at which the scale factor vanishes corresponds to tg = 
implying that the collapse ends in a finite proper time. 
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III. GRAVITATIONAL COLLAPSE WITH TACHYON AND BAROTROPIC FLUID 



The model we employ to study the gravitational collapse has an interior space-time geometry as described above 
(see e.g. (TK^UIIHI)- Concerning the matter content, we consider a spherically symmetric homogeneous tachyon field 
together with a barotropic fiuid. We use an inverse square potential for the tachyon field given by [29H3T| . 



F(0) = 



(3.1) 



where Vq is a constant. Note that we can consider the potential (3.11 as two (mirror) branches (see figure [T]), upon 
the symmetry —0, treating the (f) > and 4> < separately. 

The total energy density of the collapsing system is therefore p — + p)^, with 



3H^ 



1^6 



(3.2) 



where pb is the energy density of the barotropic matter, whose pressure in terms of the barotropic parameter, 7, 
satisfies the relation Pb = (7 — l)Pb, the barotropic parameter being positive, 7 > 0. Then, the Raychadhuri equation 
for the collapsing system can be written as 



-2H = 



7Pb- 



(3.3) 
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Furthermore, the equation of motion for the tachyon field can be obtained as 



The energy conservation of the barotropic matter is 



Ya 

V 



Pb + 37iJpb = 0. 



(3.4) 



(3.5) 



Integrating (3.5), the energy density for barotropic matter can be presented as pb = Pob^ '^''j where pob is a constant 



of integration. The total pressure is given by 



P = P0 + (7 - l)Poba 



(3.6) 



In order to have a physically reasonable matter content for the collapsing cloud, the tachyon field and the barotropic 
fluid would have to satisfy the weak and dominant energy conditions^ [T|[51|3Sj: for any time- like vector y*, the energy- 
momentum tensor Tik satisfies TikV^V^ > and T^'^Vk must be non-space-like. Since for the symmetric metric (2.1) 
we have T/ — —p{t) and — Tg — — p{t), then, for the energy densities and the pressures of the tachyon field 
and the barotropic fiuid the energy conditions amount to the following relations 



P0 > 0, 
Ph > 0, 



P4,+P4,> 0, 
Ph+Ph> 0, 



\Ph 



I < P0> 
< Ph- 



(3.7) 
(3.8) 



The first two expressions in equations (3.7) and (3.8 1 denote the 'weak energy condition' (WEC), and the last 
expression denotes the 'dominant energy condition' (DEC). 

For a tachyon field, it is straightforward to show that the WEC is satisfied: 



> 0, and p0 + p0 



> 0. 



(3.9) 



In addition, since 0^ < 1, then V{(I)){1 — cjP') ^/"^ > V{4))\J 1 — 02 from which it follows that \p^\ < p^, and hence, the 
DEC is satisfied for the tachyon matter. 

For the barotropic fluid considered herein, we assume a regular initial condition for the collapsing system, such as 
P06 > 0, indicating that pf, = pofc'J^^^ > 0- In addition, since Pb = (7 — l)pb, then Pb + Pb = IPb > provided that 
7 > 0; through this paper we will consider the positive barotropic parameter 7. Hence, the WEC is satisfied. On the 
other hand, considering the DEC {pb > \pb\) it follows that the barotropic parameter 7 must satisfy the range 7 < 2 
(cf . the phase space analysis in section III A ) . 



A. Phase space analysis 



We analyze in this (sub)section the gravitational collapse of the above setup by employing a dynamical system 
perspective. We introduce a new time variable t (instead of the proper time t present in the conioving coordinate 
system {t,r,9,ip}), defined as 

r^-lo,[^)\ (3.10) 

where < r < 00 (the limit t — corresponds to an initial condition of the collapsing system (a — >■ oq) and the limit 
r — >■ cxD corresponds to a — > 0). For any time dependent function / = f{t), we can therefore write* 



dr 3H' 



(3.11) 



^ We note that satisfying the weak energy condition implies that the null energy condition (NEC) is held as well. 
* We recall that throughout this paper H < 0, i.e., d < is assumed. 
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We further introduce a new set of dynamical variables x, y and s as follows 

V _ pb 



y 



(3.12) 



From the new variables (3.11), (3.12) and the system (3.3)-(3.5) with (3.1), an autonomous system of equations is 
retrieved: 



dx 
dr 
dy 
dr 
ds 
d^ 



A 



= /i 

= /2 

= h = s{l - s){j - x') 



Note that A is given by 



A: 



V 



1/3/2 



(3.13) 
(3.14) 
(3.15) 

(3.16) 



which for the pot enti al (3.1 ) brings A — ±2/vn) as a constant; for the > branch, A > 0, and for < branch, 
A < 0. Equation (3.2) in terms of the new variables reads^. 



y 



s = 1. 



(3.17) 



The dynamical variables have the range y > 0, s < 1, and — 1 < x < 1. This brings an effective two-dimensional 
phase space^. 

Setting (/i, /2, /3)|(£Ce,i/c,Sc) = '^an obtain the critical points {xc,yc,Sc) for the autonomous system. The 

stability can be subsequently discussed by using the eigenvalues of the Jacobi matrix A, defined at each fixed point 

{•^c: '^c^7 as 



/ 9fi dfi dh\ 
I Fix dy da \ 



A 



OX dy as 

dh dh dl 

dx dy ds 

1 dfs dj3 dfs 

\ dx dy ds 



(3.18) 



' \{Xa,Vc 



Solutions, in terms of the dynamical variables, in the neighborhood of a critical point g"'* can be extracted by making 
use of 



q^{t) 



(3.19) 



with the perturbation Sqi given by 



Sq^ = ^(<Zo)iexp(OiV), 



(3.20) 



where qi = {x,y,s}, Q are the eigenvalues of the Jacobi matrix, and the {qo)l are constants of integration. We have 
summarized the fixed points for the autonomous system and their more relevant stability properties in table |l] 



Table|T]is complemented wi th fig ure [2| The autonomous system equations ( 3.13 )-( 3.15) unfolds a three dimensional 
phase space constrained by (3.17). Therefore the main fixed points (and the trajectories nearby) are located on 
surfaces of this three dimensional phase space. Furthermore, from table |Tj we can focus our attention on three 
situations characterized by having s = 0, s = 1 and s = sq. Briefly, in more detail: 



^ Or equivalently, the surface where the trajectories will be present can be written as j/^ = (1 — s)'^ (l — x-^). 
® In the absence of a barotropic fluid, the effective phase space is one-dimensional. 



7 



point 


X 


y 


s 


Existence 


Stability 


(a), [a] 


1 








for all A; 7 < 1 
for all A; 7 > 1 


Stable node 
Saddle point 


(b), [b] 


-1 








for all A; 7 < 1 
for all A; 7 > 1 


Stable node 
Saddle point 


(c), [c] 




yo 





for all A; 7 > 7i 
for all A; 7 < 7i 


Unstable node 
Saddle point 


(d), [d] 








1 


for all A; 7 / 
for all A; 7 = 


Saddle point 
Unstable node 


(e) 




37 


So 


A > 0, 7 < 71 < 1 


Unstable node 


[e] 




37 


So 


A < 0, 7 < 71 < 1 


Unstable node 


(/), [/] 


1 





1 


for all A; 7 > 1 
for all A; 7 < 1 


Stable node 
Saddle point 


(5), [g] 


-1 





1 


for all A; 7 > 1 
for all A; 7 < 1 


Stable node 
Saddle point 


(h), [h] 


1 





Sl 


7 = 1 


Stable node 




-1 





Si 


7 = 1 


Stable node 



Table I: Summary of critical points and their properties for the </> < (denoted with square brackets, i.e., as, for example, [a]) and > 
branches (denoted with regular brackets, i.e., as, for example, (a)). 




Figure 2: Trajectories in phase space and critical points: (i) Left plot represents the phase space region {x,y,s) constrained by | |3.17[ l. 
Therein we also depicted all the fixed points (see table[l]l, except point (e), for Vb = 4/9, 7 = 0.5. (ii) In the right plot we considered the 
conditions as Vq = 4/9, 7 = 1.5. We can illustrate that going from 7 < 1 to 7 > 1 (from left to the right plot) reverses the direction of the 
trajectories, i.e., in the left plot the vector field is directed towards points (a) or (fe). In the right plot the vector field is directed towards 
points(/) or (g). 

• Point (a): The eigenvalues are Ci — —2, (,2 = —1 and Cs = (7 ^ !)■ For 7 < li all the characteristic values are 
real and negative, then the trajectories in the neighborhood of this point are attracted towards it. Hence, (a) 
is a stable node (attractor). Finally, for 7 > 1, all characteristic values are real, but one is positive and two are 
negative, the trajectories approach this point on a surface and diverge along a curve: this is a saddle point. 

• Point (b): The eigenvalues are Ci = —2, (2 = ^1 and — il ^ !)• This point has the same eigenvalues of 
point (a) and similar asymptotic behavior, being also a stable node for 7 < 1 (see the left plot in figure [s]) and 
a saddle point for 7 > 1. 

• Point (c): This fixed point has eigenvalues (i = 0, C2 = i^^Uo + 6j/o)/6 > and Cs — il ~ ^^Vo/'^): which 
all are real and yo = — A^/6 + -^/l + For 7 > 71 = two components are positive. However, 
from a numerical investigation we can assert that this corresponds to an unstable saddle (see left plot in figure 
[3]). On the other hand, for 7 < 71, one component is negative and other is positive, and hence, a saddle point 
configuration would emerge. 
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• Point (d): The eigenvalues are = 1, C2 = —7 and Cs = —7- As 7 > 0, trajectories approach this point on a 
surface (the in-set) and diverge along a curve (the out-set). This is a saddle point (see left plot in figure [2]). 

• Point (e): This point is located at (— -^7, §3-, sq), where so = (1 ^ 

get 



< 1 



37 



< 1 



A2yi^)- Fi'oni the constraint < sq < 1 we 

(3.21) 



I.e., 



< 7 < 71 < 1 



(3.22) 



When 7 — )■ 0, points (e) and (d) become coincident (see table When we consider sq = 0, we obtain 7 = 71, 
X = ~^\/yo and y = z/q. Consequently point (e) and (c) become coincident. Therefore, point (e) can be 
found along a curve that joins points (c) (on the surface with s — 0) and (d) (on the surface with s = 1). The 
eigenvalues are Ci — 0, with ^2 and C3 given by, 



C2,3 =^ (2 - 7 ± v/(l-7) (4-16so7)+72) 



(3.23) 



The eigenvalues (3.23) are non negative for 7 < 71 and for 7 = 71, C2 > and C3 — 0. 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 3: (i) The top left plot represents a section, of the three dimensional phase space presented in figure |2] labeled with s = 0. The 
dashed and dotted lines are the zeros of eqs. l |3.13| l- |"3.15| l. The full line represents constraint | |3.17| . Therefore, the fixed points are found 
on their intersections. We considered the conditions as Vq = 4/9, 7 = 0.5. In this example we can locate the position of fixed points (a), 
(6) and (c) according to table[l] (ii) The center plot represents a section, showing fixed point (e), labeled with s = 0.76. The dashed (and 
also dotted) lines are the zeros of eqs. | |3.13[ l-| pi.l5[ l and the full line represents constraint | |3.17[ l, as before. We considered the conditions 
as Vb = 4/9, 7 = 0.5. (iii) Finally, in the top right plot we represent a section, labeled with s = 1. We can identify fixed points (g), (d) 
and (/). We considered the conditions as Vb = 4/9, 7 = 1.2. 



In the plot at the center of figure [3] we provide a section of the phase space showing point (e), where it can be 
seen that all trajectories are divergent from it. This behavior is characteristic of an unstable node and expected 
whenever point (e) is found along the curve that connects (c) and (d). 

• Point (/): The eigenvalues are Ci = —2, (,2 = —7 and Cs = (1 ^ 7)- For 7 > 1, all the characteristic values are 
real and negative, then the trajectories in the neighborhood of this point are attracted towards it. Hence, (/) 
is a stable node (attractor). Finally, for 7 < 1, all characteristic values are real, but one is positive and two are 
negative. The trajectories approach this point on a surface and diverge along a curve; this is a saddle point. 

• Point (g): The eigenvalues are Ci — —2, (2 — ~7 and Ca = (1 ~ 7)- This point has the same eigenvalues of 
point (/) and similar asymptotic behavior, being also a stable node for 7 > 1 (see the right plot in figure |3]) and 
a saddle point for 7 < 1. 

• Point (h): This point is found along the line segment with x = 1, y = and s £ ]0, 1[. The eigenvalues are 
Ci = —2, C2 — —1 and Ca — 2(si — 1), where si G ]0, 1[. In this case, all the characteristic values are real and 
negative, then the trajectories in the neighborhood of this points are attracted towards it. Hence, {h) is a stable 
node (attractor). 
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• Point (i): This point is found along the hne segment with x = —1, y = and s € ]0, 1[. The eigenvalues are 
Ci = —2, C2 = — 1 and ^3 = 2(si — 1). This point has the same eigenvalues of point (h) and similar asymptotic 
behavior, being also a stable node. 

Before proceeding, let us add two comments. On the one hand, note the transition that occurs, when going from 7 < 1 



to 7 > 1, an intermediate state as discussed (by means of a numerical study) in subsection IIIB3 This bifurcation 
behaviour is made explicit through the numerical methods employed. These allowed us to confirm the results on 
the dynamical system analysis and to further assess in regions like the transition from tachyon dominance to fluid 
dominance, verifying all possible scenarios with those two tools. On the other hand, figure |2] deserves some attention 
when we consider the trajectories approaching the stable nodes (a) — [h) (left plot) and (/) — [g) (right plot). Therein, 
the trajectories approach a; — ±1 along a segment line containing the stable nodes. This situation implies that the 



energy density (3.2) diverges after reaching a point where x — ±1, y — > and s — > si. This observation will be 



important in the next sections concerning the possible outcomes of the gravitational collapse. 



B. Analytical and numerical results 



Let us herewith discuss this section possible outcomes regarding the collapsing system, employing elements from 
both our analytical as well of numerical study. 



1. Tachyon dominated solutions 
From the trajectories in the vicinity of (a) and (6), attractor solutions can be described. 

The asymptotic behavior of s near the point (a) can be approximated as s ~ Sc + exp(— r) = exp(— r); hence, as 
T — >■ cx) (i.e., a — >■ 0), s vanishes. Moreover, the time derivative of the tachyon field is given by ~ 1, that is, the 
tachyon field (j){t) has a linear time dependence and can be approximated as (see figure |4]) 

(j)(t) ~ i + (3.24) 




/ t t 

Figure 4: Behavior of the area radius, tachyon field and its time derivative over time. We considered the initial conditions as: t^ = 0, 
POb = li = 4/9, a(0) = ao, (0) = —0.6 and (f>o = 1. We also have 7 = 0.1 (full lines) and 7 = 1.2 (dashed lines). 



It should be noted that, for the 4> > branch, 0o is positive on the initial configuration of the collapsing system 
(where t = 0) and hence, the tachyon field increases with time, proceeding downhill the potential. Within a finite 
amount of time, the tachyon field reaches its maximum^ but finite value (f>{ts) = (ps at tg < (j>oi with the minimum 
(but non-zero value) V oc t^J^. As the tachyon potential decreases, the dynamical variable y = 7^ vanishes. 

Thus, as (/) — > 1, the energy densityof the system diverges. Furthermore, the tachyon pressure = —V{(f>){l~(p)^ 
vanishes asymptotically^ (see figure for plots of the energy density, barotropic pressure and tachyon field pressure) . 



^ The time at which the collapse reaches the singularity is finite. Thus the tachyon field at the singularity remains finite as 0(is) = ts+(f>o 

® For the (f> < branch, 0o is negative at the initial condition. Thus, the absolute value of the tachyon field starts to decrease from the 
initial configuration as <j>{t) = t — \(f)o\ until the singular point at time ts < \<t>o\, where tachyon field reaches its minimum but non-zero 
value 4>s. This leads it uphill the potential until the singular epoch, where the potential becomes maximum but finite. 
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0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 

t t t 



Figure 5: The energy density, barotropic pressure and tachyon field pressure. We considered ti = 0, a(0) = ao, Vb = 4/9 and 0o = 0.6 
witii 7 < 1 (full lines) and 7 > 1 (dashed lines). The total effective pressure (p(, + p^) is divergent and negative (for 7 < 1) in the final 
stage of the collapse. 




Figure 6: Behavior of the ratios ^^^3 (left) and (right). We considered the initial conditions as: ti = 0, po6 = li = 4/9, a(0) = aoi 
0(0) = —0.6 and t/iQ = 1- We illustrate the transition between tachyon and fluid dominated solutions as a function of the barotropic 
parameter. In these collapsing cases the potential does not diverge and the tachyon field does not vanish. 



The total energy density of the cohapsing system is given approximately by the energy density of tachyon field (see 
figure 6 for numerical solutions), which corresponds to a dust-like matter near the singularity as oc 1/a'^ . From 
this relation for the tachyon energy, we may write the energy density of the system near point (a), 

p « pofl-^ (3.25) 

This induces a black hole at the collapse final state. 

As far as (6) is concerned, the time dependence of the tachyon field can be obtained as 

0(t) ~ -t + (/)o. (3.26) 

For the </> > branch, the tachyon field decreases from its initial value at i = (where (/)o > 0), moving uphill the 
potential. Then, the potential will reach a maximum (but finite) value when the tachyon field reaches its minimum 
and nonzero value (i.e. 0s as t —> is < (/)o). In the limit case ts — (t>Q, the tachyon field vanishes and the potential 
diverges. I.e., y and s vanish, the Hubble rate increases faster than the potential and the barotropic fluid energy 
density diverges, that is, the total energy density of the system asymptotically diverges.^. Likewise, when — > 1 and 
a — 0, the tachyon matter behaves as dust matter. The fate of the collapse for this fixed point is as well a black hole 
formation. 

The asymptotic solutions provided by the fixed points (a) and (6) correspond to a dust-like solution with a vanishing 
pressure for the tachyon field, whose energy density reads oc 1/a^. This is consistent with the WEC and DEC for 
the tachyon matter being satisfied, as was mentioned before. Concerning the status of the energy conditions for the 
barotropic fluid, as we indicated before, regularity of the initial data for the collapsing matter respects the WEC. On 
the other hand, stability of the solution in this case ensures that 7 < 1, which satisfies the sufficient condition for the 
DEC. 



On the other hand, for the < branch, the tachyon field increases from its initial condition as time evolves, proceeding to ever less 
negative values as </< — > 0~. In this case, it proceeds downhill the tachyon potential till the system reaches a = at t = ts, where the 
Hubble rate and hence the total energy density of the system diverge. This implies that the time at which the collapse system reaches 
the singularity is always ts < </>o- 
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2. Fluid dominated solutions 

From the trajectories in the vicinity of (/) and {g), solutions can be described with some having an attractor 
behaviour. 

The asymptotic behavior of s near the point (/) can be approximated as s w 1 + exp(— r); hence, as r — > oo (i.e., 
a — > 0), s 1. Moreover, the time derivative of the tachyon field is given by ~ 1; that is, the tachyon field 4>{t) has 
a linear time dependence and can be approximated as (see figure |4]), 

0(0 ~ i + (3.27) 

The total energy density of the collapsing system is given approximately by the energy density of the barotropic 
fluid (see figure p| where pb/ (3i?^) 1 while p^j,/ (3H^) — > for 7 > 1 ), which goes, near the singularity, as 
Ph cx From this relation for the fluid energy, we may write the energy density of the system near point (/), 

p « pobfl^'^. (3.28) 

For 7 < 2/3, the ratio F/R — ir^pofeO.^"^''', converges as the singularity is reached leading to the avoidance of trapped 
surfaces, but since the corresponding flxed point (/) turns to be a saddle, then the resulting naked singularity is not 
stable. For 2/3 < 7 < 1 the ratio F/R diverges and the trapped surfaces do form. But still the point (/) is saddle 
and the resulting black hole is not stable. The case 7 > 1 corresponds to a stable solution for which the ratio F/R 
goes to infinity as the collapse advances. Then, the trapped surface formation in the collapse takes place before the 
singularity formation and thus the final outcome is a black hole. 

As far as (g) is concerned, the time dependence of the tachyon field can be obtained as 

(l){t)~-t + (f,Q. (3.29) 

For the </> > branch, the tachyon field decreases from its initial value a.t t — (where ipQ > 0), moving uphill the 
potential, the potential will reach a maximum (but finite) value when the tachyon field reaches its minimum and 
nonzero value (i.e. — ^ (/)s as t — )• < 0o)- In the limit case tg = 4>o, the tachyon field vanishes and the potential 
diverges. I.e., y and s vanish, the Hubble rate increases faster than the potential and the barotropic fluid energy density 



diverges, that is, the total energy density of the system, is given by equation (3.28), and asymptotically diverges. 



Similar to point (/) the mass function of the system for the fixed point solution {g) is given by F/R = ^r^ poi,a^~^'^ , 
which diverges for 7 > 1. Therefore, the resulting singularity in this case will be covered by a black hole horizon. 

As far as the energy conditions are concerned for the fixed point solutions (/) and (g), we find that the tachyon 
fleld satisfies the WEC. Also the DEC remains valid case as well. On the other hand, for the barotropic fiuid, the 
WEC is satisfied initially and will hold until the endstate of the collapse. The stable solution in this case corresponds 
to the range 7 > 1 which satisfles DEC as well. 



3. Tracking solution: black hole and naked singularity formation 



In this subsection, we will discuss a difi^erent type of solutions, where the fiuid and tachyon appear with a tracking 
behaviour. Let us introduce this situation as follows. 

Interesting and physically reasonable tracking solutions can be found, where cf) — >■ ±1, when we consider 7 — >■ 1, 
i.e., a situation whereby the emergence of points (h) and (i) will be of relevance as attractors. The transition from 
tachyon dominated to fluid dominated scenarios, like those described in previous sections, is not straightforward. In 
this situation, the tachyon field and barotropic fiuid compete to establish the dominance in the late stage of the 
collapse. In figure [6] we have a illustration of this kind of solutions. The mentioned dominance seems to depend 
strongly on the initial ratio p,p/pb at an earlier stage of the collapse and also on the value of 7. Such a dependence on 
the initial conditions can lead to a set of solutions between those provided by fixed points (a), (b) (tachyon dominated 
solutions) and by points (/), (g) (fluid dominated solutions). From a dynamical system point of view, this corresponds 
to have trajectories asymptotically approaching x — > ±1 in sections where s is between and 1 (see figure [2]). At 
the end of the collapse we observe that 3^ ^ 3^ < 1, as illustrated in figure |6] Therefore, in this scenario, the 
trajectories would convey a collapsing case in which the energy density of the tachyon field and of the barotropic fluid 
are given by p^ (x ph oc a~^^ . This shows a tracking behavior for the collapsing system [501133]. Moreover, the total 
energy density of the collapse, in terms of a, reads 



p (X a-3T. 



(3.30) 
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Equation (3.30) shows that the energy densities of the tachyon field, the barotropic matter and hence, their total 
(for the collapsing system), diverges as a — > 0. The ratio of the total mass function over the area radius is given by 



F 
R 



1 



oc 



(3.31) 



Equation (3.31) subsequently implies that, for an adequate choice of values {'^,4>q), trapped surfaces can form as the 
collapse evolves and a few scenarios can be extracted. More precisely, for the range 7 > |, for both cf) > and (j) < 
branches, the final fate of the collapse is a black hole. For the case in which 7 < |, the ratio F/R remains finite as the 
collapse proceeds and an apparent hor izon i s delayed or fails to form; the final state is a naked singularity (a solution 
for the choice of '— ' sign in equation (3.23 1). The tracking solution indicates 7 = | as the threshold (illustrated in 
figure [3| , which distinguishes a black hole or a naked singularity forming. 

Therefore, under suitable conditions, we can determine whether it is possible to have the formation of a naked 
singularity. In fact, if we assume a very unbalanced i nitial ratio Pip/ Pb with p^^ <^ pob ^^nd a barotropic fluid having 
7 < |, then we can have a situation where the ratio (3.31) is converging. The set of initial conditions described by 
Po</> ^ Pob are equivalent to consider the barotropic fluid as initially dominant. If this specific unbalanced distribution 
of matter is allowed to evolve into a regime where the tachyon dominates then the system will evolve until becomes 
comparable to pi,. The singularity is reached in finite time and it can happen before the tachyon can effectively 
dominate. In figure [t] we have a graphical representation of the ratio F/R. It can be seen that the ratio F/R remains 
finite for 7 < 2/3, while the energy density is diverging, as the collapse proceeds and apparent horizon is delayed or 
fails to form till the singularity formation. As the right plot shows, the validity of WEC is guaranteed throughout 
the collapse scenario for both barotropic fluid and tachyon field. Also the DEC is valid for the solutions that exhibit 
naked singularity, i.e., those for which 7 < 2/3. 




Figure 7: In the left plot it is shown the ratio of the total mass function over the area radius. We considered Vb = 4/9 and po6 = 1100. 
It is shown that for 7 > 2/3 we have black formation and for 7 < 2/3 we have naked singularity formation. The right plot shows the weak 
energy condition over time for the barotropic fluid (b.f) and for the tachyon field (0). Therein we considered ti = 0, a(0) = ao, Vb = 4/9, 
^ = 1100 and = 0.6 with 7 = 0.6. 



4- Other solutions 

Point (c) does not represent an attractor and hence its vicinity will not corresponds to a collapse endstate. However, 
trajectories emerging from it would proceed towards (a) or (h). The time dependence of the tachyon field near (c) 
can be obtained integrating = as 




(3.32) 



In the (/) > branch, the tachyon field starts to increase with time from its initial value at f = (where (t){t = 0) = 
(/)o > 0), proceeding uphill the potential. Then, a stage is obtained where a maximum (but finite) value is reached, 
when the tachyon field goes to its minimum and nonzero value at < \^/3Jyo. Asymptotically, y — — > yo and 

s ~ 7^ — > 0. Since > is a constant, the Hubble rate at t — reads H = —y^^VjSyo and the total energy density 
of the system is essentially the energy density of the tachyon field and is given by 

P^3H2 (3.33) 

yo 
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which is proportional to the tachyon potential. Therefore, at the time ts < X^/sJyo, when the tachyon field reaches its 
non-zero minimum, the total energy density of the system remains finite. In the limit case, when the time ts = A-^/s/yo, 
the tachyon field vanishes, the tachyon potential diverges. Alternatively, we can also write that tachyon behaves as 

« a—. (3.34) 
When a — ^ 0, the tachyon field vanishes, and hence, the potential and energy density of tachyon field in equation 



(3.33), as a function of a, is given by p oc V{4>) a ^ ^° and diverges at a = 0. The dynamics of the system in 
vicinity of this fixed point is given by a; w "^a/^O' U ^ Vo + (o-o/o.)^'^^ and s w (ao/a)"^^^, where ^2 > 0. 

Concerning (d), the energy density of the tachyon field is approximately given by the tachyon potential ~ ^(0) 
(with = 0) and thus, in its vicinity, the tachyon potential remains constant and finite. Furthermore, y — V/3H^ ~ 
and pb/iH^ — 1 imply that the total energy density is given by the energy density of barotropic matter only: 
p ~ 3H^ ~ pb- This solution diverges as a — 0. 

The fixed point (e) for the range 7 < 71 presents solutions for the collapsing system in which asymptotically is 
given hy 4> — — ^7, for the cj) > and 4> — \/l for the </) < branch. For the > branch, the time dependence of 
the tachyon field is given by 

0W = -\/7i + </'o. (3.35) 

The tachyon starts its evolution from the initial configuration 0o > 0, and decreases uphill the potential. Then, at 
the time t ~ tg = ^ the tachyon field vanishes, the potential diverging^°. 

A thorough numerical analysis of the trajectories in the vicinity of point (e) has shown that it behaves as an unstable 
node. We can observe in the center plot of figure [3] an illustration of the general behavior near this point, namely 
having trajectories diverging from it. 

C. Exterior geometry 

In order to complete the full space-time geometry for the herein collapsing model, we need to match the homogeneous 
interior space-time to a suitable (inhomogeneous) exterior geometry. For a perfect fluid gravitational collapse set up, 
with equation of state p = (7 — 1) p, the pressure does not necessarily vanish at the boundary. E.g., matching the 
internal geometry filled with matter (and radiation) to a boundary layer (which is crossed with the radiation), which 
could in turn be matched to an exterior geometry, not completely empty (e.g., filled by radiation). More concretely, 
matching the interior with a generalized exterior Vaidya space-time across the boundary given by r = r;, |28j . 

Similarly to [34], we proceed by considering the interior metric which describes the collapsing cloud as 

ds^ = -dt^ -f a^{t){dr^ ^ r^dn^) (3.36) 

and the exterior one in advanced null coordinates {v, R) given by 

dsl = -f{v, R)dv^ + 2dvdR + R^dn^. (3.37) 

In these coordinates the outermost boundary of trapped surfaces is simply given by the relation f{v, R) = 0. It 
should be noticed that in general the formation of the singularity at a = is independent of matching the interior 
to the exterior space-time. In order to find a suitable exterior metric function f{v,R) we resort to an Hamiltonian 
perspective of the model. The total Hamiltonian constraint is given by |35j 



H{a, TTa, TT^) ^ Ha + Htf, + 7ih 

yV^2 + ^-6^2 + po,a-3(7-i), (3.38) 



a 



12a 

where tTq and tt^ are the conjugate momentums for the scale factor a and for the tachyon field 0, respectively. 



Furthermore, the Hamilton equations for parameters a, tTq and tt^ can be obtained by using the equation (3.38) as 



For the < branch, the time dependence of the tachyon field is given by <f>{t) = ^y^t + ipQ. Therefore the tachyon field starts its 
evolution from the initial configuration at rj>o < 0, decreasing in time, going uphill the potential. 
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follows: 



an 



m 

da 

cm 

' d<j) 



^+3(7-l)poba-(i+3(--i» 



12o2 



+ Til 



(3.39) 
(3.40) 
(3.41) 



On the other hand, since the Hamiltonian constraint, "H = 0, must be held across the boundary E, equation (3.38) 
reduces to 



"^^144^+^°*^" +^°^ 6^ 



(3.42) 



whereby, substituting for tTq and tt^ from equation (3.39)- (3.41), we get 



1 



(3.43) 



Furthermore, it is required for the junction condition rba(t) — R{t) and the equation of motion for scale factor r^d = R 
to be satisfied at the boundary of two regions. By substituting these conditions in the Hamiltonian constraint equation, 
^{{VjR) — 0, given by (3.43), we get 



(1-/) 



9(1 - /) + 6p, 



'Ob 



4-37' 



Plh 



-67 



(3.44) 



Notice that, we have substituted the four- velocity of the boundary being seen from the exterior by i? = — -^1 — / (cf. 
[34|). Then, by solving equation (3.44 1 for /, the boundary function is obtained simply as 



fiR) = 1 - 



V 



1-6^ 



1 



rf,\ 1+3(7-1) 



(3.45) 



Let us now study the behavior of the boundary function for the stable fixed point solutions we obtained in the 
previous subsections. We note that the above expression is valid when both the tachyon field and barotropic fluid 
are present, thus in the regimes where the tachyon field is dominant the contribution of the fluid to Hamiltonian 
constraint is set aside. For the dust-like solutions, descri bed by points (a) and (fe), the energy density of tachyon field 
at the bounda ry r;,, included in second term of equation (3.45 ), is given by ~ Porf/R^. So, the boundary function 
/ in equation (3.451 reduces to: 



fiR) = 1 - 



M 



where 



47r , 



(3.46) 



and Po = Po/27rr^ = const. Equation (3.46) for f{R) constrains the exterior space-time to have a Schwarzschild metric 
in advanced null coordinates, providing an interpretation of the collapsing system as a dust ball with the radius rj, 
and the density po- On the other hand, for the fixed point solution (e), by using the equation (3.301 in (3.451, we 
have 



f{R) = 1 - 



M 



where 



M=-\po 



2so + l 
So 



Pob ^ = const. 



(3.47) 



where 37 — 2 > (notice that, this corresponds to the case in which 7 = 71 > 2/3 in previous subsection). 

In order to get a possible class of dynamical exterior solutions, we proceed by considering the following metric at 
the boundary [M] 



ds' = -[l- 



2M{R, v) 
R 



dv^ + 2dvdR + r^dn^ 



(3.48) 
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where 

M(i?,«)=rn(.)-^^^-^|^, (3.49) 

is the total mass within the collapsing cloud. Matching the for exterior metric function gives 

2m(w) g{v) M 



1 



R (27-3)i?27-2 i?37-2 ' 

2rfi{v) g{v) 



(3.50) 



R (27 - 3)i?2'r-2 ' 

in which, the second part stands for matching for the extrinsic curvature. Differentiation of the first expression in 



(3.50) with respect to time and using the second one, we get 

2m{v) 2{j-l)g{v) M{3-/-2) 



(3.51) 



i?2 (27-3)ii27-l 3i?37-l ■ 

Multiplying the first expression in equation (3.501 by R^^ and after adding the result with the above equation, we 
get 

g(.)^- ^^;;^^ (3.52) 



Now, by substituting for g{v) into equation (3.50), we obtain 



M7 

= 6(3-27)fl(.)37-3 ■ (3.53) 



Then, as seen from equation ( 3.31[ ) for 7 < 2/3, where the trapping of light has failed to occur, the exterior geometry 



is dynamical in contrast to the tachyon dominated regime in which a space-like singularity forms with a static exterior 
space-time. Finally, for fluid dominated solutions, depending on the value of 7, both naked singularities and black 
holes may form, the mass being different to those of tracking solutions and the exterior geometry being static or 
dynamical, respectively. 

IV. CONCLUSIONS, DISCUSSION AND OUTLOOK 

In this paper we considered a particular setting among the spherically symmetric class of models for gravitational 
collapse, with a tachyon field and a barotropic fluid as matter content. We restricted ourselves to the marginally bound 
case (cf. [TJ ?ITi ,36j for a description and details) . The tachyon potential was assumed to be of an inverse square form. 
Our objective was to establish (i) which final states would occur (namely, a black hole or a naked singularity), (ii) how 
each matter component will compete (the fluid being conventional, whereas the tachyon bringing some workable but 
intrinsic non-standard effects from string theory) and, (iii) which will eventually be the determinant component at the 
end of the collapse process. More precisely, in our opinion it is of interest to assert, if and how, at later stages, effects 
induced by the tachyon (a scalar field (among others) found in string theory context), could allow interesting features 
to b e ev entually discussed (i.e., what asymptotic behaviour emerges). As far as the tachyon field is concerned, from 



eq. (3.4) we can have that with —1 < (f) < 1 and H <0, (j) terms act like an anti-friction contribution at a collapsing 



phase, within an uphill evolution for the (j) 2 potential for tachyon field (i.e., when < 0). Moreover, other terms 
(i.e. the term including ^) would have an anti-friction effect as well. 

Determining therefore the outcome of the gravitational collapse in our system^^, i.e., whether, e.g., a black hole 
or naked singularity would form, was not a straightforward assessment. We considered an analytical description by 
means of a phase space analysis |201 132j. discussing several asymptotic behaviours; these were also subject to a careful 
study involving a numerical investigation, which added a clearer description of the possible dynamical evolutions. 



The case of a standard scalar field was investigated in |21l |36] 
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Within this setting, for a spatiahy homogeneous interior space-time, we found a situation where the tachyon was 
dominant, with a black hole forming^^. A cosmological framework involving only a FRW geometry, with a tachyon 
field and a barotropic fluid, was investigated in [50], focusing on the late time (dark energy like) stages. It is interesting 
to note that while for a tachyon dominated regime, an inflationary-like (i.e., accelerated expansion) scenario leads 
to violation of the strong energy condition (SEC)^"^ [20 , in our collapse scenario we have instead a corresponding 
asymptotic stage, where tachyon dominance leads to black hole formation, satisfying the SEC at the flnal state of the 
collapse (cf. figure Q). The same behaviour is observed in the fluid dominated regime, where for 7 > 1 the SEC 
is satisfled by the fluid, with a black hole formation. However, for the tracking solutions, where a naked singularity 
forms, i.e., pob >> P041 and 7 < 2/3, the SEC is violated. 

Moreover, in further comparison with the set of fixed points found in }20j for an expanding (accelerating) universe 
{H > 0), we have found analytically two additional critical points (/ and g) for the collapse process which correspond 
to a barotropic dominated collapsing regime which ends in a black hole for 7 > 1, respecting the WEC and SEC. 
Indeed, for the late-time acceleration of the universe, filled with a tachyon field and a fluid, the phase space analysis 
in [20] predicted two class of solutions: a tachyon dominated solution; and a tracking solution [20]. Rather differ- 
ently, corresponding solutions herein our paper turn to be unstable towards the singularity in the collapse process. 
Nevertheless, two other solutions are also provided in our collapsing scenario: a tachyon dominated solution where 7 
satisfles the range 7 < 1, initially; and a barotropic dominated solution for which the barotropic parameter holds the 
range 7 > 1. All these solutions predict a black hole formation as collapse end state. 

Further regarding the tracking solutions (of a cosmological nature) indicated in |20j . In our collapsing system, a 
different and rather interesting set of states is found, within the context of tracking behavior for the barotropic fluid 
plus tachyon field content. Being more concrete, these solutions have that a black hole or naked singularity forms. In 
particular, in this situation, we found that it is possible to define by a numerical appraisal in particular, the threshold 
7 = 2/3, separating black hole and naked singularity solutions for the gravitational collapse. We have also discussed 
the specific conditions leading to the formation of a naked singularity. Therefore, we concluded, that in our model, if 
the collapse starts with an unbalanced distribution of the matter content favoring the barotropic ffuid, i.e. <C p;,, 
then towards the final stage of the collapse they evolve until ~ ph, when a naked singularity forms. However, since 
the NEC is satisfied, the strong curvature condition along the null geodesies can be preserved and the singularity can 
be strong in the sense of [S7] . 

We think it is fair to indicate that we employed a potential for the tachyon, whereas for </> — ^ 0, the tachyon 
should not induce a divergent behavior as far as string theory advises [38H42I . In fact, an exponential- like potential 
for the tachyon could bring a richer set of possible outcomes [221 |3S1 221 HI] , including a better behaved and possibly 
a more realistic evolution when dealing with — 0. 

Finally, let us add that it will be of interest to investigate (i) other scenarios for the geometry of the interior 
space-time region, (ii) specific couplings between the tachyon and the ffuid, within e.g., a chamaleonic scenario for 
gravitational collapse and black hole production (broadening the scope in [iSHTf] ). (iii) adding either axionic, dilatonic 
or other terms (e.g., curvature invariants) that could be considered from a string setup, but at the price of making 
the framework severely less workable and (iv) whether explicit quantum effects can alter the outcomes presented in 
this paper. To this purpose latter, we plan to use ingredients brought from loop quantum gravity (cf. ref. ^48j ) . 

V. ACKNOWLEDGMENTS 

The authors would like to thank P. Joshi, A. Khaleghi, C. Kiefer, F. C. Mena, A. Vikman, and J. Ward for making 
useful suggestions and comments. They are also grateful to the referee for the useful comments and suggestions on 
the issue of energy conditions. YT is supported by the Portuguese Agency Fundagao para a Ciencia e Tecnologia 
through SFRH/BD/43709/2008. This research work was also supported by the grant CERN/FP/109351/2009 and 
CERN/FP/116373/2010. 



[1] P. Joshi, Gravitational Collapse and Space-Time Singularities, (Cambridge University Press, 2007). 



It is worthwhile to compare the result obtained herein for a homogeneous (tachyonic) collapsing matter field with the gravitational 
collapse of a k-essence scalar field with non-standard kinematic terms in [30) . where the scalar field has a dependence to the radius r 
i.e., (p = <^(r). In both models, the collapsing systems lead to the black hole formation. 
Satisfying the SEC for the tachyon field demands that + 3p^ > 0. Thus, the SEC holds if 0^ > 2/3. 



17 



[2] S. W. Hawking, G. F. R. Ellis, The Large Scale Structure of Space- Time, (Cambridge University Press, 1974). 
[3] R. Penrose, Riv. Nuovo Cimento 1 252 (1969). 
[4] T. Harada, Phys. Rev. D 58 (1998) 104015; 

T. Harada and H. Maeda, Phys. Rev. D 63 (2001) 084022; 

R. Goswami and P. S. Joshi, Class. Quantum Grav 19 (2002) 5229; 

R. Giambo, F. Giannoni, G. Magli and P. Piccione, Gen. Rel. Grav 36 (2004) 1279; 

J. F. Villas da Rocha and A. Wang, Class. Quantum Grav 17 (2000) 2589; 

R. Giambo, F. Giannoni, G. Magli and P. Piccione, Class. Quantum Grav 20 (2003) 4943. 
[5] P. Szekeres and V. Iyer, Phys. Rev. D 47 (1993) 4362; 

S. Barve, T. P. Singh, and L. Witten, Gen. Rel. Grav. 32 (2000) 697; 

A. A. Coley and B. O. J. Tupper, Phys. Rev. D 29 2701 (1984); 

K. Lake, Phys. Rev. D 26 (1982) 518. 
[6] R. Giambo, Class. Quantum Grav 22 (2005) 2295; 

S. Bhattacharya, R. Goswami and P. S. Joshi, Int. J. Mod. Phys. D 20 (2011) 1123; 

S. Bhattacharya, arXiv:1107.4112 [gr-qc]. 
[7] A. Ori and T. Piran, Phys. Rev. Lett 59 (1987) 2137; 

A. Ori and T. Piran, Phys. Rev. D 42 (1990) 1068; 

T. Foghzzo and R. Henriksen, Phys. Rev. D 48 (1993) 4645. 
[8] F. C. Mena, B. C. Nolan and R. Tavakol Phys. Rev. D 70 (2004) 084030; 

S. G. Gosh and N. Dadhich, Gen. Rel. Grav. 35 (2003) 359; 

T. Harko, S. K. Cheng, Phys. Lett. A 226 (2000) 249. 
[9] C. Gundlach and J. M. Martm-Garcfa, Living Rev. Relativity 10 (2007), 5. 

F. E. Schunck and E. W. Mielke, Class. Quant. Grav. 20 (2003) R301; arXiv:0801.0307 [astro-ph]. 
[10] R. Giambo, F. Giannoni and G. MagU, J. Math. Phys.49 (2008) 042504; Gen. Rel. Grav. 41 (2009) 21. 
[11] K. Ganguly and N. Banerjee, Gen. Rel. Grav. 43 (2011) 2141. 
[12] A. I. Janis, E. T. Newman and J. Winicour, Phys. Rev. Lett 20 (1968) 878; 

M. Wyman, Phys. Rev. D 24 (1981) 839; 

B. C. Xanthop oulos and T. Zannias, Phys. Rev. D 40 (1989) 2564; 
O. Bergmann and R. Leipuik, Phys. Rev. 107 1157 (1957); 

H. A. Buchdahl, Phys. Rev 115 (1959) 1325. 
[13] M. Choptuik, Phys. Rev. Lett. 70 (1993) 9. 
[14] D. Christodoulou, Ann. Math. 140 (1994) 607. 

[15] F. I. M. Pereira and R. Chan, Int. J. Mod. Phys. D 17 (2008) 2143; arXiv:0803.2628 [gr-qc]. 
[16] M. Jankiewicz and A. A. Sen, arXiv:0602.085 [gr-qc]. 

[17] A. Mazumdar, S. Panda and A. Perez-Lorenzana, Nucl. Phys. B 614 101 (2001). 

[18] Yun-Song Piao, Rong-Gen Cai, Xinmin Zhang, and Yuan-Zhong Zhang, Phys. Rev. D 66 (2002) 121301; 

M. Fairbairn and M. H. G. Tytgat, Phys. Lett. B 546 1 (2002). 
[19] A. Sen, Int. J. Mod. Phys. A, 20 (2005) 5513; arXiv: 0410.103 [hep-th]. 
[20] Juan M. Aguirregabiria and Ruth Lazkoz, Phys. Rev. D 69, 123502 (2004). 
[21] R. Goswami, P. Joshi and P. Singh, Phys. Rev. Lett. 96, 031302 (2006). 
[22] S. A Hayward, Phys. Rev. D 53, 1938 (1996); 

C. W. Misner and D. H. Sharp, Phys. Rev. 136, B571 (1964); 

D. Bak and S. J. Rey, Class. Quantum Grav. 17, L83 (2000). 
[23] T. P. Singh, Phys. Rev. D 58, 024004 (1998). 

[24] E. Malec, N. 6 Murchadha , Phys. Rev. D 47, 1454 (1993). 
[25] R. Goswami, and P. Joshi, [arXiv: 0410.144 [gr-qc] (2004). 

[26] A. H. Ziaie, K. Atazadeh and Y. Tavakoh, Class. Quant. Grav. 27, 075016 (2010); arXiv: 1003.1725 [gr-qc]; 

A. H. Ziaie, K. Atazadeh, S. M. M. RasouU (2011), arXiv: 1106.5638 [gr-qcj. 
[27] P. S. Joshi and R. Goswami, Class. Quantum Grav. 24 2917 (2007); 

R. Goswami, P. S. Joshi, C. Vaz and L. Witten, Phys.Rev. D 70 084038 (2004); 

P. S. Joshi and R. Goswami, arXiv:0711.0426vl [gr-qc] 
[28] P. S. Joshi and I. H. Dwivcdi, Class. Quant. Grav, 16 41 (1999); A. Wang and Y. Wu, Gen. Relat. Grav, 31 107 (1999). 
[29] Mohammad R. Garousi, Nucl. Phys. B 584, 284 (2000); E. A. Bergshoeff, M. de Roo, T.C. de Wit, E. Eyras, S. Panda, 

J. High Energy Phys. 05, 009 (2000). 
[30] R. Akhoury, D. Garfinklea, and R. Saotomea, J. high Energy Phys. 1104, 096 (2011); arXiv: 1103.0290 [gr-qc]. 
[31] Andrei Frolov, Lev Kofman, Alexei Starobinsky, Phys. Lett. B 545 8 (2002); arXiv: 0204187 [hep-th]. 
[32] H. K. Khalil, Nonlinear Systems, 2nd edition (Englewood Cliffs. NJ: Prentice Hall, 1996); Jack Carr, Applications of Center 

Manifold Theorem (Springer- Verlag, 1981); J. Guckenheimer, P. Holmes, Nonlinear Oscillations, Dynamical Systems and 

Bifurcation of Vector Fields (Springer- Verlag, 1983). 
[33] A. R. Liddle, and R. J. Scherrer, Phys. Rev. D 59, 023509 (1999). 
[34] B. K. Tippett and V. Hussain, Phys. Rev. D 84 (2011) 104031. 

[35] A. A. Sen Phys. Rev. D 74 043501 (2006); M. Tsamparlis and A PaUathanasisa, arXiv:1111.5567 [gr-qc]. 
[36] R. Goswami, and P. Joshi, arXiv: 0410.144 [gr-qc]. 



18 



[37] F. J. Tiplcr, Phys. Lett. A 67 (1977) 8; F. J. Tiplcr, C. J. S. Clarke, and G. F. R. Ellis, in General Relativity and 
Gravitation, edited by A. Held (Plenum, New York, 1980), vol. 2, p. 97; C. J. S. Clarke, The Analysis of Space-time 
Singularities (Cambridge University Press, Cambridge, England, 1993). 

[38] A. Feinstein, Phys. Rev. D 66, 063511 (2002). 

[39] Z-K. Guo, Y. S. Piao, R. G. Cai, and Y. Z. Zhang, Phys. Rev. D 68, 043508 (2003). 
[40] L. R. W. Ambramo and F. Finelli, Phys. Lett. B 575, 165 (2003). 

[41] E. J. Copeland, M. R. Garousi, M. Sami, and S. Tsujikawa, Phys. Rev. D 71, 043003 (2005). 

[42] I. Quiros, T. Gonzalez, D. Gonzalez, Y. Napoles, R. Garci'ar-Salcedo, and C. Moreno, Class. Quant. Grav. 27, 215021 

(2010). 

[43] A. Sen, .]. High Energy phys. 12, 021 (1998). 

[44] A. A. Sen, J. high Energy Phys. 04, 048 (2002); A. A. Sen, J. high Energy Phys. 07, 065 (2002); 
16, 281 (1989); 

S. H. Ghate, R. V. Saraykar, and K. D. Patil, Pramana, J. Phys. 53, 253 (1999). 

[45] V. Folomccv, arXiv: 1205.2974 [astro-ph]. 

[46] V. Folomeev and D. Singleton, Phys. Rev. D 85 (2012) 064045; arXiv:1112.1786 [astro-phj. 
[47] V. Folomeev, Phys. Rev. D 85 (2012) 024008; arXiv: 1108.3395 [astro-ph]. 

[48] Yascr Tavakoli, Joao Marto, Amir Hadi Ziaie and Paulo Vargas Moniz, Semiclassical gravitational collapse with tachyon 
field and barotropic fluid, to appear in PRD (2013). 



